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We investigate, in an exact manner, the phase structure of the micromaser system in terms of the physical 
parameters at hand like the atom cavity transit time r, the atom-photon frequency detuning Aui, the number of 
thermal photons rib and the probability a for a pump atom to be in its excited state. Phase diagrams are mapped 
out for various values of the physical parameters. At sufficiently large values of the detuning we find a "twinkling" 
mode of the micromaser system. A correlation length is used to study fluctuations close to the various phase 
transitions. 
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There are few systems in physics which exhibit a rich 
struture of phase transitions that can be investigated un- 
der clean experimental conditions and which at the same 
time can be studied by exact theoretical methods. The 
micromaser system, which is a remarkable experimental 
realization of the idealized system of a two-level atom 
interacting with a second quantized single-mode electro- 
magnetic field, provides us with such an example (for 
reviews and references see e.g. [|J). 

Many features of the micromaser are of general inter- 
est. It can e.g. be argued that the micromaser system 
is a simple illustration of the conjectured topological ori- 
gin of second-order phase transitions 12] . Various aspects 
of stochastic resonance can, furthermore, be studied in 
this system M . The micromaser also illustrates a feature 
of non- linear dynamical systems: turning on randomness 
may led to an increased signal to noise ratio B. 

In a typical realization of the micromaser, the pump 
atoms which enter the cavity are assumed to be prepared 
in an incoherent mixture, i.e. the density matrix pa of 
the atoms is diagonal with diagonal matrix elements a 
and b such that a + b = 1 . In terms of the dimensionless 
atomic flux parameter N — R/j, where R is the rate 
injected atoms and 7 is the damping rate of the cavity, 
the stationary photon number probability distribution is 
then well known H and is given by 



Pn = Po 



n.a 



rib m + Naq n 



■ (1 + n b )m + Nbq n 



Here q m = q(m/N) and 
q{x) = 



x + A 2 



(eVx + a 2 ) 



(i) 



(2) 
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where we have defined the dimensionless detuning A = 
Aio / (2g\/N) and pump 8 = grVN parameters. Fur- 
thermore, g is the single photon Rabi frequency at zero 
detuning of the Jaynes-Cummings ( JC) model |J . The 
overall constant po is determined by Y^Lo P' n = 1 ■ 

The theory as developed in Refs. ||[7| suggests the ex- 
istence of various phase transitions in the large N limit as 
the parameter 9 is increased. A natural order parameter 
is then the average photon "density" (x), where ( ) de- 
notes an average with respect to the distribution Eq. (|l|) 
and x = n/N . An exact treatment, in the large N limit, 
of the micromaser phases structure and the correspond- 
ing critical fluctuations in terms of a conventional cor- 
relation length has been given in Refs. ||. Spontaneous 
jumps in (n)/N and large correlation lengths close to mi- 
cromaser phase transitions have actually been observed 
experimentally [n|M| . Most of the theoretical and exper- 
imental studies have, however, been limited to the case 
a = 1 and A = 0. It is the purpose of the present paper 
to study the phase structure of the micromaser system 
for general a and A using methods which are exact in the 
large N limit ||. As will be argued below, several new 
intriguing physical properties of the micromaser system 
are then unfolded. 

In order to obtain an appropriate expression for p n 
which describes the various micromaser phases, we notice 
that the equilibrium distribution in Eq. (|l|) can be re- 
written by using the Poisson summation technique p0[ . 
The equilibrium distribution then takes the form 



^>=W^r A ' V( " 



(3) 



where we have defined an effective potential V(x) — 
Er=-co y fc(a;) and 

V k (x) =- J dv ln[tu(i/)] cos{2TiNkv) , (4) 

Jo 



w(x) = 



n b x + aq(x) 
(1 + n b )x + bq(x) 



(5) 



In the large N limit Eq. (g|) can be simplified by mak- 
ing use of a saddle-point approximation. The saddle- 
points are then determined by Vq(x) = 0. If such non- 
trivial saddle-points exist we say that they describe maser 
phases. In a topological analysis of the second-order 
phase transitions || of the micromaser system, Vq (x) will 
play the role of a Morse function. 

If the only global minimum of Vq (x) corresponds to x — 
0, we expand the effective potential around the origin. 
The micromaser is then in a thermal phase and we obtain 



Pa 



n b 



'eff 



l + n b + b9l ff 



(6) 



which is normalizable provided 8\i * (a — 6) < 1, where 
we have defined B\ ff = sin 2 (0A)/A 2 . If a < (1 + A 2 )/2 
the distribution Eq. (g) is always valid. If A = nor- 
malizability requires that 8 2 (a — b) < 1, i.e. 9 must be 
sufficiently small if a > 1/2. If A ^ the mean value of 
photons obtained from Eq. (|6|) will be a periodic function 
of A0 with a maximum (n b + a/A 2 )/(l + (b- a) /A 2 ) for 
A0 = (n + l/2)7r, where n = 0, 1, ... . The corresponding 
minimum is nj, and occurs for A9 — nix. In the thermal 
phase the micromaser can therefore be in a "twinkling" 
mode for A / 0, i.e. the mean number of photons can 
exhibit periodic variations as a function of 9. 

Non-trivial saddle-points of the effective potential 
Vo(x), which exist only if a > (1 + A 2 )/2, can be para- 
metrically represented in the form B 



x + A 2 = (a - b) sin 2 cf> 
9 l ' ; 



\Ja — b | sin <j>\ 



(7) 



with 



> 



= arcsin(|A|/Va — b). These saddle-points 
do not depend on n b . If, for a given 9, there are several 
saddle-points, the actual maser phase is described by the 
saddle-point which corresponds to the global minimum 
ofVbO). 

The first extremum of Vo(x) occurs at 0q = 9{<j>o). 
Critical points, where new extrema of Vq(x) appear, are 
determined by Vq(x) = 0, i.e. non-trivial solutions of 
tan 4> = (f>. Corresponding to the positive solutions 
<f> = (f>k, where k = 1,2, ... , we have the critical pump 
parameters 9k = 4>k/{\ sin0fc|Va — b). 

Using the substitution <j) = 9\J x + A 2 , the potential 
Vq(x) in Eq. (0) can now be written in the form 



A| 



Vo{4>,6) 

fxj) ln[ 



n b + aq(4>,9) 



l + n b + bq(4>,9) 



(8) 



where q(4>, 9) = 9 sin 4>/<f> ■ By choosing cj) correspond- 
ing to Eq. (fn) the effective potential is always at an ex- 
tremum. We then get a multi-branched effective poten- 
tial Vq — Vq(6). Branch k is swept out when (f> varies 
in the range <fio + kit < <j) < (k + l)n — (f>Q, where 



k = 0, 1, 2, ... . Except for the first branch (k = 0), each 
of these branches is doubled- valued. One sub-branch 
then corresponds to a maximum (Vq(x) < 0), which is 
swept out first as </> increases, and the other corresponds 
to a minimum (Vq'(x) > 0). For a given branch k, these 
sub-branches coincide at the critical point 9k (see e.g. 
Fig. 0). 
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FIG. 1. The extremal values of the effective potential Vo(0) 
when a = 1, Ub — 0.15, A = (upper figure) and |A| = 0.5 
(lower figure). The critical parameters 9q and 8k are shown 
as well as the pump parameters 91^+1 f° r co-existence of two 
maser phases. 8l k+1 , marked by a circle, and 9* k correspond 
to first-order transitions. 

In Fig. |TJ we illustrate such branches as well as the cor- 
responding value of the order parameter (a;) . The tran- 
sition from the thermal phase to the first maser phase at 
8 = 0g is always a second-order phase transition. The 
transition from one maser branch k to the neighboring 
maser branch k + 1, both with Vq(x) > 0, corresponds 
to a first-order phase transition and occurs at 9 = 9% k+1 . 
For various k such critical points are marked by circles 



in Fig. [j]. As A increases the nature of these transitions 
changes. In Fig. [j] with |A| = 0.5 the first two maser 
phases e.g. never intersect. Instead they intersect with 
the thermal phase. This phenomena corresponds to the 
twinkling behavior in the thermal phase even though we 
do not have a strictly periodic behavior. The twinkling 
phenomena will now, however, be more pronounced since, 
in the large N limit, the maser will be "dark" in the ther- 
mal phase. For a = 1 and rib = 0.15 the phase separation 
of the first two maser phases occurs at |A| « 0.408. With 
increasing value of A 2 < 1 more phases are separated. 

Critical transition lines can now be determined by con- 
sidering intersections between the various micromaser 
phases. Hence, we can construct a phase diagram in e.g. 
the a-0-parameter space. The result is shown in Fig. g 
for the same parameters as in Fig. [j]. In the upper phase 
diagram in Fig. |2| the various micromaser phases are well 
separated, i.e. they are uniquely separated by the criti- 
cal lines. The first critical line is in general determined 
analytically by the condition B\^ [a — b) = 1. The other 
critical lines are determined by means of analytical and 
numerical methods. In the lower phase diagram of Fig. 
with |A| = 0.5, micromaser phases which appear for suf- 
ficiently large values of 8 will again be well separated. 

The order parameter {n)/N can now also be studied as 
a function of e.g. the probability a. Fig. y shows such a 
plot. When we cross transition lines in the phase diagram 
the order parameter will make discrete jumps and hence 
also exhibits a plateau-like behavior. 

Let us now consider long-time correlations in the large 
TV limit as was first introduced in Refs. S. These 
correlations are most conveniently obtained by making 
use of the continuous time formulation of the micro- 
maser system [ |ll| . The vector p formed by the diag- 
onal density matrix elements of the photon field then 
obeys the differential equation dp/dt = —jLp, where 
L = Lc — N(M — 1). Here Lc describes the damping of 
the cavity , i.e. 



(L C )r 



(rib + l)[n5„, r 
■n b [{n + l)8 n . 



(n + l)<5„+i, m ] 
■ nS niTn+1 } , (9) 



and M = M(+) + M(-) , where M(+) nm = 
bq n +iSn+i, m + a(l - q n +i)S n ,m and M(-) nm = 
aq n Sn,m+i + b(l — q n )5 n ,m have their origin in the JC- 
model |||J. The lowest eigenvalue Ao = of L then 
determines the stationary equilibrium solution p — p° as 
given by Eq. (|lj). The next non-zero eigenvalue A of X, 
which we determine numerically, will then determine 
typical scales for the approach to the stationary situa- 
tion. The joint probability for observing two atoms, with 
a time-delay t between them, can now be used in order 
to define a correlation length "f A (t) ||. At large times 
t — > oo, we define the atomic beam correlation length £4 

by[pl 
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FIG. 2. Phase diagrams for the same set of parameters as 
in Fig. hi The upper figure is discussed in the main text. 
The dashed curves in the lower figure correspond to the lines 
determined by d^ff ( a ~b) = 1- When not visible they overlap 
with the solid critical lines. The first critical line corresponds 
to a second-order (thermal-maser) transition. For any other 
critical line the transition is first-order (second- order) to the 
left (right) of its minimum (determined by sin 2 (A(9) = 1) 



unless it intersects with another critical line, 
are indicated by circles. 



Triple points 



(10) 



which then is determined by A, i.e. 7^4 = 1/A. For 
photons we define a similar correlation length £@. It 
follows that the correlation lengths are identical, i.e. 
(,A = (,c = £, H- The correlation length 7^ is shown in 
Fig. [| for I A| = 0.5. In the large N limit, the clear peaks 
in the correlation length occur at the critical pump pa- 
rameters 6*o, 9* k and 0* kk+1 . In the thermal phase we then 
have 7£ ~ l/(l + (o— 6)6* 2 ^), provided we are not to close 
to a critical line. When we move from the thermal phase 
to a maser phase the correlation length will increase. At 
the critical line (a — b)9^ ff = 1, the correlation length 

will behave as ("f£,)crit — (a — b)y/N/(a + (a — b)rib)/2, 
apart from a weak dependence of the detuning parameter 

A. 
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the various phase transitions and critical fluctuations 
of the micromaser system at non-zero detuning and for 
pump atoms prepared in a statistical mixture. We have 
revealed new novel features of the system as e.g. the 
plateaus in the number of photons (n)/N as a function 
of a as well as a new twinkling mode of the micromaser 
system for large detuning. This twinkling phenomena 
has a close resemblance with the observed revivals of the 
micromaser system W . 
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FIG. 3. Plateaus of the order parameter (a;) = {n)/N ver- 
sus a when A — 0, nb = 0.15 and 9 = 25 for various values of 
N. Each steg_ corresponds to a point on one of the transition 
lines in Fig. 



At the transition lines 6 = 9 



fefe+i 



the methods of |H| can 



be used to show that the correlation length depends ex- 
ponentially on the combination N{a + n;,(a — b)). One 
can then also show that, as A 2 approaches its maximal 
value a — b for saddle-points to exist, 7^ — > 1 exponen- 
tially fast. In view of the recent experiments on trapping 
states in the stationary field configuration of the micro- 
maser p2| , we notice that high peaks occur in the atomic 
correlation length for such states ||. The dependence of 
these trapping state peaks on the physical parameters at 
hand can be studied using the methods presented here. 



20 



15 



10 



ex 



0: 



e: 



9* 







0= 1 


- log{~iO t 


N = 125 


n b = 0.15 

|A| = 0.5 




N = 25 




o 5 


10 


15 20 



e 



FIG. 4. The logarithm of the correlation length 7^ as a 
function of 8 for various values of N(25, 50, ...125). The ver- 
tical lines indicate the critical values of 6 (see Fig. hi) . 

In conclusion we have studied, in the large N limit, 
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